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Abstract. A linear system of differential equations describing a joint motion of 
thermoelastic porous body and incompressible thermofluid occupying porous space 
is considered. Although the problem is linear, it is very hard to tackle due to the 
fact that its main differential equations involve non-smooth oscillatory coefficients, 
both big and small, under the differentiation operators. The rigorous justification 
is fulfilled for homogenization procedures as the dimensionless size of the pores 
tends to zero, while the porous body is geometrically periodic. As the results, we 
derive Biot's like system of equations of thermo-poroelasticity, system of equations 
of thermo-viscoelasticity, or system of non-isotropic Lame's equations depending 
on ratios between physical parameters and geometry of porous space. The proofs 
are based on Nguetseng's two-scale convergence method of homogenization in pe- 
riodic structures. 
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Introduction 

In the present publication we consider a problem of a joint motion of 
thermoelastic deformable solid (thermoelastic skeleton), perforated by a sys- 
tem of channels (pores) and incompressible thermofluid occupying a porous 
space. We refer to this model as to model (NA). In dimensionless variables 
(without primes) 

x' = Lx, t' = rt, w' = Lw, 9' = 

TV* 

the differential equations of the model in a domain fiel 3 for the dimen- 
sionless displacement vector w of the continuum medium and the dimension- 
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less temperature 9, have a form: 

a Tp-Q^ = div^P + pF, (0.1) 

99 d 
a ^v~Ql = div ^^^) ~ a e —div x w + (0.2) 

P = xa^fx, ^) + (1 - x)«a©(x, «;)-(? + tt)I, (0.3) 

ff = p + ^ + *M, (0.4) 

p + x a pdiv x w = 0, (0.5) 
7r + (1 - x)(a v div x w - a 9s 0) = 0. (0.6) 

Here and further we use notations 

B(x,u) = (l/2) (V x u + (V x u) T ), 

P = XPf + 0-- X)Ps, Cp = xcpf + (1 - x)c ps , 

a* = xa><f + (1 - x)®xs, a e = xa ef + (1 - x)®e s - 

In this model the characteristic function of the porous space x{x) is a known 
function. For derivation of (jO.ip - fl0.6p and description of dimensionless 
constants (all these constant are positive ) see [TT] . 

We endow model (NA) with initial and boundary conditions 

dw 

w\ t=0 = w , -77-\t=o = v , 9\ t=0 = 9 , xett (0.7) 
at 

w = 0, 9 = 0, xeS = dtt, t>0. (0.8) 

From the purely mathematical point of view, the corresponding initial-boundary 
value problem for model (NA) is well-posed in the sense that it has a unique 
solution belonging to a suitable functional space on any finite temporal inter- 
val (see [H]). However, in view of possible applications this model is ineffec- 
tive. Therefore arises a question of finding an effective approximate models. 
If the model involves the small parameter e, the most natural approach to 
this problem is to derive models that would describe limiting regimes arising 
as £ tends to zero. Such an approximation significantly simplifies the original 
problem and at the same time preserves all of its main features. In the model 
under consideration we define e as the characteristic size of pores I divided 
by the characteristic size L of the entire porous body: 

I 

£ = T 
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But even this approach is too hard to work out, and some additional 
simplifying assumptions are necessary. In terms of geometrical properties 
of the medium, the most appropriate is to simplify the problem postulating 
that the porous structure is periodic. Further by model (NB) £ we will 
call model NA supplemented by this periodicity condition. Thus, our main 
goal now is a derivation of all possible homogenized equations in the model 
(NB) £ . 

We accept the following constraints 

Assumption 0.1. domain Q = (0, l) 3 is a periodic repetition of an elemen- 
tary cell Y £ = eY, where Y = (0, l) 3 and quantity \je is integer, so that Q 
always contains an integer number of elementary cells Yf. Let Y s be a "solid 
part" of Y, and the "liquid part" Yf - is its open complement. We denote 
as 7 = dYf fl dY s and 7 is C 1 -surface. A porous space is the periodic 
repetition of the elementary cell eYf, and solid skeleton Q £ s is the periodic 
repetition of the elementary cell eY s . A boundary T £ = dfl £ s fl dfl'j is the 
periodic repetition in fl of the boundary £7. The "solid skeleton" Q s is a 
connected domain. 

In these assumption 

X(x) = x £ (x) =x(x/e) , 

c p = c £ p (x) = x £ (x)c pf + (1 - x £ (x))c ps , 

p = p £ ( x ) = x £ (x)pf + (1 - X £ (x))p s , 

a M = a%{x) = x £ {x)a ><f + (1 - x £ {x))a„ s , 

a e = a £ e {x) = x £ {x)uej + (1 - x £ {x))a es , 

where x{v) i s a characteristic function of Yf in Y. 

We say that a porous space is disconnected (isolated pores) if 

7 n BY = 0. 

In the present work we suppose that all dimensionless parameters depend on 
the small parameter e and there exist limits (finite or infinite) 

UmaJe) = fi , lima\(e) = A , limo; r (e) = r , limQ;„(e) = p*. 

e\0 e\0 e\0 e\0 

Moreover, we restrict ourself with the case when r < 00 and 
/i — 0, = 00, < A < 00. 
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If t = oo,then, re-normalizing the displacement vector and temperature 
by setting 

w — > a T w, 9 — > a T 6 

we reduce the problem to the previous case. The condition = oo means 
that the liquid in consideration is incompressible. 

Using Nguetseng's two-scale convergence method [9j [12] we derive Biot's 
like systems of thermo-poroelasticity or system of non-isotropic Lame's, de- 
pending on the ratios between dimensionless parameters and geometry of the 
porous space. 

Different isothermic models have been considered in [TJ], [2], [13], [HEBE], 

§1. Model (NB) £ 

§1. Formulation of the main results. 

As usual, equations (I0.1l) - (l0.6p are understood in the sense of distribu- 
tions. They involve the equations (10.11) - (10.61) in a usual sense in the domains 
fl £ j: and Q £ s and the boundary conditions 

[#] = o, [w] = 0, x g r e , t > 0, (1.1) 

[P] = 0, [aiy x 9] = 0, x G T £ , t > (1.2) 

on the boundary T £ , where 

[(p](x ) = <P( B )(x ) - (p(j)(x ), 
<ff s )(x ) = lim (p(x), ipm(x ) = lim tp(x). 

X — * XQ X — > XQ 

x E 0| I 6 fij 

There are various equivalent in the sense of distributions forms of rep- 
resentation of equations (I0.1I) - (I0.2I) and boundary conditions (ll.ll) - (ll.2l) . In 
what follows, it is convenient to write them in the form of the integral equal- 
ities. 

Definition 1.1. Five functions (w £ , 8 £ , p £ , q £ , tt £ ) are called a generalized 
solution of model (NB) £ if they satisfy the regularity conditions in the 
domain Qt = x (0, T) 

w £ , B(x, w £ ), dw x w £ , q £ , p £ , 7T e , r, V x 6 £ G L 2 {Q T ) (1.3) 
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in the domain fl T = Q x (0,T), boundary conditions (10.81) . equations 

q £ = p £ + -% + x £ <y ef 6 £ , (1.4) 

p e + x e a p div x i« e = 0, (1.5) 
tt £ + (1 - ^(e^div^ - a 0s 6 E ) = (1.6) 

a.e. in Qt, and integral identities 

jf (a T pW • ^ - X e %©(^, ™ £ ) : D(x, ^) - p £ F ■ if + 
{(1 - x £ )a A 0(x, ™ E ) - (<f + vr e )I} : D(x, <p))dxdt + 



«rP £ (^o • "^l*=o ~ • V\t=o)dx = (1.7) 



for all smooth vector-functions ip = ip{x,t) such that <p\an = <p\t=T = 
d<p/dt\ t= T = and 

J ([a T c £ p 6 £ + c^div^)^ - a%V x 9 E ■ V s £ + V£)dxdt 

+ / (a r c£0§ + a £ e dw x w £ )£\ t=0 )dx = (1.8) 

for all smooth functions £ = £(«c,t) such that £|an = £|*=t = 0. 

In (II. 41) by A : B we denote the convolution (or, equivalently, the inner 
tensor product) of two second-rank tensors along the both indexes, i.e., A : 
B = tr(B*oA) = Zl =1 A 3 B Jl . 

Suppose additionally that there exist limits (finite or infinite) 



lima„(e) = vq, lim a v (e) = r] , lima a<8 (e) = x 0s , \ima ef (e) = /3 0f , 

e\0 e\0 e\0 e\0 

lim a gs (e) = (3 0s , lim ^ = fi u lim — = X/. 



e\0 "" v ' ' e\0 S 2 ' " e\0 a 



In what follows we suppose to be held 



Assumption 1.1. l)Dimensionless parameters in the model (NB) £ satisfy 
to next restrictions 

/i = 0; 0<t o + Hu x 0s , X/, A , Vo', 
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To, Xf, XOs, vq Pof, Pos A < oo. 

2) Sequences {y/a^( l -X £ )Vwp}, {y/c£v e Q }, {y/a^6 e }, {^/a^(l-x £ )div x w £ }, 
{v^X e div x «>g}, {Voac^(1-X 6 )V«§}, {yo~a~(l-x e )div x <^}, {^/d^o^diy, 
{ag}, {bg} are uniformly in e bounded in L 2 (Q) and \F\, \dF/dt\,^f, d^/dt e 
L 2 (fi T ). 

Here 

a^div^ + pF^O), 

c£bg = div a .(a&V a .0g) - a £ e div x v £ + *(aj, 0), 

¥ £ = tfa^ix, v £ ) + (1 - if)a x B(x, «)§) + 

(^(opdiva-twg + c^div^) + (1 - x e )o^div a .ti;g)I. 

In what follows all parameters may take all permitted values. For example, 
if To = or tJq 1 = 0, then all terms in final equations containing these 
parameters disappear. The following Theorems ll.ltil.2l are the main results 
of the paper. 

Theorem 1.1. For all e > on the arbitrary time interval [0, T] there exists 
a unique generalized solution of model (NB) 6 and 

maxJ||™ £ (t)|,Va-x £ ^^ (1.9) 

\\e £ hn T + V^\\x E V x e e h,n T + ||(1 - X £ )V x 6 e \\ 2 ,n T < Co, (1.10) 

h e h,n T + b e ||2,n T + — H-^rlkoT + lk e lkn r < Co (1.11) 
where Co does not depend on the small parameter e. 

Theorem 1.2. Functions w e and 8 s admit an extension u e and d e respec- 
tively from Q e sT = Q e s x (0, T) into Qt such that the sequences {u 6 } and {i) £ } 
converge strongly in L 2 (Qt) and weakly in L 2 ((0, T); W / 2 1 (f2)) to the functions 
u and d respectively. At the same time, sequences {w 6 }, {0 s }, {p 6 }, {q E }, 
and {n e } converge weakly in L 2 (Q T ) to w, 6, p, q, and it, respectively. 
The following assertions for these limiting functions hold true: 
(I) If H\ = oo then w = u, 9 = i? and the weak limits u, p, q, and n 
satisfy in Qt the initial-boundary value problem 

d 2 u 

ToP W + V(q + n)-pF= } (112) 

div x {\ A s : D(x, u) + B s {div x u - ^tf) + B{q}, 



(roe P + ^(l-m))^-^ + (al--)(^) n = div x (B e .W) + *, (1.13) 
Vo ot r] dt r] dt 

-(tt + (g) n ) + C S : B{x, u) + a s (div x u - ^(0 - (#>„)) + a?(g - (g) n ) = 0, 

(1.14) 

-(vr + (g)n) + div x u + (1 ~ m)Afe (t? - (0> n ) = 0, (1.15) 
?-(g>n=P + /Vn(0-<0) n ), (L16) 

u>/iere 

p = mp/ + (1 - m)p s , c p = mc pf + (1 - m)c ps , m = J x(v)dy. 

The symmetric strictly positively defined constant fourth-rank tensor Aq, 
constant matrices Cq,Bq Bf, strictly positively defined constant matrix B 
and constants a s Q , a\ and a\ are defined below by Eqs. (I4.32p - (I4.34j) and 
d43U). 

Differential equations (I1.12l) - (ll.l6p are endowed with initial conditions at 
t = and x E VI 

{t + (3 0s ){$-%) = 0, t (u - u ) = t (— - v ) = 0; (1.17) 

and boundary conditions 

ti(x,t) = 0, u(x,t) = 0, x e S, t>0. (1.18) 

(II) If the porous space is disconnected, then w = u and strong and 
weak limits u, d, p, q, it together with a weak limit 9* of the sequence {x £ @ 6 } 
satisfy in equations (11.121) . (11.141) - (11.151) . the state equation 

q-(q)n = P + Po f (0 f -(9%), (1.19) 

and heat equation 



+ (nc r , + ^)(1 - ,n) Tt - >^ Tt + W - -)<^> 

div x (B 9 ■ W) + *. (1.20) 

Here 9* is defined below by formulas (14.39j) - (14.44p . 

The problem is endowed with initial and boundary conditions (11. 17ft - 
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(Ill) If Hi < oo then strong and weak limits u, w$ , 0* , p, q and n of 
the sequences {u £ }, {$ e }, {x £ w £ }, {x e 9 £ }, {p e }> {q 6 } an d i 716 } satisfy the 
initial- boundary value problem in Q?, consisting of the balance of momentum 
equation 

QPujf d^u 
cHv x {\qAq : D(x, u) + Bydiv x u + Bfq}, 

where Aq ; Bq B{ are the same as in Eq. (11.121) . continuity equation (11.141) . 
continuity equation 

-(tt + (q)n) + div x w f + (1 ~ m)/ ^ _ (tf) n ) = ( m _ l)div x u, (1.22) 
Vo Vo 

state equation fll. 19|) . heat equation (11.201) and Darcy's law in the form 



dw? du 



m- 



I 1 d 2 u 
J Bi{Hi,t-T)-{-V x q + p f F-T p f —){x,T)dT (1.23) 



dt dt 
if tq > /ii > ; Darcy's law in the form 
dwf du 



dt dt 

if T o — and, finally, Darcy's law in the form 

Q w f :>»■ ' ''' 



+ B 2 (Hi)-(-V x q + p f F) (1.24) 



r i?3 -l^ + M-S 3 )- / (-V x q + p f F)(x,T)dr (1. 

at <9t T pf J 



25) 



«/>i = 0. 

T/ie problem is supplemented by boundary and initial conditions (II. 17ft - 
(I1.18P for the displacement u and temperature $ of the rigid component and 
by the boundary condition 

w f (x, t) ■ n(x) = 0, (x, t) e S = dtt, t>0 (1.26) 

for the displacement of the liquid component. In Eqs. (II. 23p - (11.261) n(x) 
is the unit normal vector to S at a point x G S, and matrices Bi(p,i,t), 
B 2 (hi)> and B 3 are given below by formulas (14. 50p - 04.551) . 

52. Preliminaries 
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2.1. Two-scale convergence. Justification of Theorems II. 1H1 .21 relies 
on systematic use of the method of two-scale convergence, which had been 
proposed by G. Nguetseng [T2] and has been applied recently to a wide range 
of homogenization problems (see, for example, the survey [9]). 

Definition 2.1. A sequence {p 6 } C L 2 (Qt) is said to be two-scale convergent 
to a limit (p G L 2 (Qt x Y) if and only if for any 1-periodic in y function 
a = o~(x, t, y) the limiting relation 

lim / p £ (x,t)a (x,t,x/e) dxdt = / / ip(x,t,y)a(x,t,y)dydxdt (2.1) 
£ ^°Jn T Jq t Jy 

holds. 

Existence and main properties of weakly convergent sequences are estab- 
lished by the following fundamental theorem pT2l [Q] : 

Theorem 2.1. (Nguetseng's theorem) 

1. Any bounded in L 2 {Q) sequence contains a subsequence, two-scale 
convergent to some limit p G L 2 (Qt x 5^)- 

2. Let sequences {p 6 } and {eV^} be uniformly bounded in L 2 (Qt)- Then 
there exist a 1-periodic in y function tp = p(x,t,y) and a subsequence {p 6 } 
such that ip, V y p G L 2 (Qt x Y), and (p £ and eV x ip e two-scale converge to ip 
and V y p, respectively. 

3. Let sequences {p £ } and {V x p £ } be bounded in L 2 (Q). Then there exist 
functions <p G L 2 (Qt) andip G L 2 (fiyxF) and a subsequence from {p £ } such 
that ip is 1-periodic in y, V y ip G L 2 (Qt x Y), and p £ and V x p £ two-scale 
converge to <p and V x p(x,t) + V ' yipi^x^t^y) , respectively. 

Corollary 2.1. Let o G L 2 (Y) and o~ £ (x) := a(x/e). Assume that a se- 
quence {p 6 } C L 2 (Qt) two-scale converges to p G L 2 (Qt x Y). Then the 
sequence a e p £ two-scale converges to op. 

2.2. An extension lemma. The typical difficulty in homogenization 
problems while passing to a limit in Model (NB) £ as e \ arises because 
of the fact that the bounds on the gradient of displacement V x w £ may be 
distinct in liquid and rigid phases. The classical approach in overcoming 
this difficulty consists of constructing of extension to the whole Q of the 
displacement field defined merely on Q s . The following lemma is valid due 
to the well-known results from [JJ, [7] . We formulate it in appropriate for us 
form: 
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Lemma 2.1. Suppose that assumptions of Sec. 1.2 on geometry of periodic 
structure hold, ip e G Wl(Vt £ s ) and ip £ = on 5| = dQ £ s D <9£7 m t/ie trace 
sense. Then there exists a function o £ G W^iP) such that its restriction on 
the sub-domain fl £ s coincide with ip £ , i.e., 

(1 - x £ (x))(a £ (x) - i/f(x)) = 0, x G Q, (2.2) 

and, moreover, the estimate 

\\o- £ h,n < C||V> e |kn S , \\^^ £ h,n < C||V x #|| 2 ,m (2.3) 

hold true, where the constant C depends only on geometry Y and does not 
depend on e. 

2.3. Friedrichs— Poincare's inequality in periodic structure. The 

following lemma was proved by L. Tartar in [TH Appendix]. It specifies 
Friedrichs-Poincare's inequality for e-periodic structure. 

Lemma 2.2. Suppose that assumptions on the geometry of Q'y hold true. 

o 

Then for any function ip *ElW\ (fij) the inequality 

I \y\ 2 dx < Ce 2 i \V x (p\ 2 dx (2.4) 

holds true with some constant C, independent of e. 

2.4. Some notation. Further we denote 

1) ($) y = J y Sdy, ($> y/ = J Yf My, ($) Ya = f Ys $dy. 

2) If a and b are two vectors then the matrix a ® b is defined by the 
formula 

(a ® b) • c = a(b • c) 

for any vector c. 

3) If B and C are two matrices, then B £§) C is a forth-rank tensor such 
that its convolution with any matrix A is defined by the formula 

(B®C) :A = B(C:A) 

4) By we denote the 3 x 3-matrix with just one non-vanishing entry, 
which is equal to one and stands in the i-th row and the j-th column. 

5) We also introduce 

J ij = + V 1 ) = ~(e* ® ej + ej ® e*), 
where (ei,e2,es) are the standard Cartesian basis vectors. 
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§3. Proof of Theorem O 

Under restriction r > estimates (11. 91) - (11. 101) follow from 

dw £ dw £ 
max J^ divj— t 2 * + V«a (t) 2) n§ 

o<t<T vv "' at B at " 

d 2 w £ dw £ d6 £ 

+ va^||-^2-(*)||2,n + v/a^lldiv^^— (t)|| 2 ,nj + vW||-~r(*)lkn) 



. — II £V7 9 2 iu e <9 2 ™ e C , . M 

+ V^llx V x -^-|| 2 ,n T + V a Ax dw x —- r \\ 2 ,n T < (3.1) 



'0_ 

where Co is independent of e. Last estimate we obtain if we differentiate 
equations for w £ and 9 £ with respect to time, multiply first equation by 
d 2 w £ /dt 2 , second equation-by d6 £ /dt, integrate by parts and sum the result. 
The same estimate guaranties the existence and uniqueness of the generalized 
solution for the model (NB) £ . 

Estimate (11.111) for pressures follows from integral identity (JT7TJ) and 
estimates (13. ip as an estimate of the corresponding functional, if we re- 
normalized pressures, such that 



/ (q £ (x,t) + n £ (x,t))dx = 0. 
Jn 



Indeed, integral identity (11.71) and estimates (13. ip imply 
(q £ + ir £ )div x il>dx\ < C|| V^|| 2 ,n- 



Choosing now if) such that (q £ + ir £ ) = div x if) 

we get the desired estimate for the sum of pressures (q £ + it 6 ). Such a 
choice is always possible (see [8]), if we put 

ip = Vip+ipo, div^Vo = 0, Aip = q £ +ir £ } ip\ m = 0, (V<p+if) )\dn = 0. 

Note that the re-normalization of the pressures (q £ +ir £ ) transforms continuity 
and state equations (jl.4j) - (jl.6p for pressures into 

q£=p£ + ^ L p % + x£{adf6£ + 1 f ) > (3 ' 2) 

— p £ + x £ div x w £ = --(3 £ x £ , (3.3) 
a p m 

—ir £ + (1 - X £ )(div x w s - ^9 £ + 7s £ ) = 0, (3.4) 
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where 



(3 £ = ((1 - X £ )div x w £ ) n , mi ) = (q £ ) n - a ef ( X £ 6 £ )n, 
(1 - = ~(q £ )n + —((!- X £ W £ )n - P £ - 

Note that the basic integral identity (jl.7p permits to bound only the sum 
(q e + n e ). But thanks to the property that the product of these two functions 
is equal to zero, it is enough to get bounds for each of these functions. The 
pressure p £ is bounded from the state equation (13.21) . if we substitute the term 
[a v I 'a p )dp s I 'dt from the continuity equation (13. 3p and use estimate (13. ip . 

Estimation of w £ and 9 £ in the case tq = is not simple, and we outline 
it in more detail. As usual, we obtain the basic estimates if we multiply 
equations for w £ by dw £ /dt, equation for 9 s by 8 s , sum the result and then 
integrate by parts all obtained terms. The only two terms F ■ dw £ /dt and 
^ • 6 £ heed additional consideration here. First of all, on the strength of 
Lemma 12.11 we construct an extension u £ of the function w £ from Q £ s into 
Q £ f such that u £ = w £ in Q £ s , u £ G W 2 l {ti) and 

C 

\\u £ hv < C\\V x u%, n < -=||(1 - X £ )V^^xW £ \\ 2 ,n- 



After that we estimate ||iu e ||2,n with the help of Friedrichs-Poincare's 
inequality in periodic structure (lemma \2.2\\ for the difference (u £ — w £ ): 



\w £ \\ 2 ,n < \\u £ \\ 2 ,n + \\u £ - w £ \\ 2 ^ < \\u £ \\ 2 , n + Ce\\x £ V x (u £ - w £ )\\ 2 , 

_ i 

< \\u £ h,n + Cs\\V x u e \\ 2> n + C{ea il 2 )\\x £ ^/^N xW £ \\ 2 ^ 

C -i 

< -=||(1 - X £ )V^V x w £ \\ 2in + C(e ail 2 )\\x £ ^V x w £ \\ 2 , n . 



The same method we apply for the temperature 9 £ : there is an extension 
§ £ of the function 9 £ from Q £ s into such that $ £ = 9 £ in Q £ s , § £ e W 2 l (Vt) 
and 

Whv < C\\VJ £ \\ 2 v < -^=||(1 - X e )v^V^1l2,n, 



> £l \ 2 ,n < -^\\(l-X £ )V^sV x 9 £ \\ 2 , n + C(eaJ)\\ X £ V^V x 9 £ \\ 2 , n . 
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Next we pass the derivative with respect to time from dw e /dt to p e F and 
bound all obtained new terms in a usual way with the help of Holder and 
Grownwall's inequalities. 

The rest of the proof is the same as for the case To > 0, if we use the 
consequence of (13. ip : 

d 2 w e 

max _— -(f) 2n < Cb. 

o<t<T at 2 

§4. Proof of Theorem 11.21 

4.1. Weak and two-scale limits of sequences of displacement, 
temperatures and pressures. On the strength of Theorem 11.11 the se- 
quences {0 £ }, {p e }, {q 6 }, {tt £ } and {w £ } are uniformly in e bounded in 
L 2 {Vtx)- Hence there exist a subsequence of small parameters {s > 0} and 
functions 9, p, q, t\ and w such that 

9 s —>0, p e — > p, q £ — >■ q, n e — >■ 7r, w £ ^ w 

weakly in L 2 (f2r) as e \ 0. 

Due to Lemma I27T1 there is a function tt e G £°°((0, T); Wj(fi)) such that 
u e = w e in Q s x (0,T), and the family {u 6 } is uniformly in e bounded 
in L°°((0, T); W 7 ^^))- Therefore it is possible to extract a subsequence of 
{e > 0} such that 

w £ -> u weakly in L 2 ((0, T); Wjftfi)) 

as £ \ 0. 

Applying again the same lemma 12.11 we conclude that there is a function 
$ e G L 2 ((0,T); ^(O)) such that ^ = fl £ in fi s x (0,T), and the family {tf £ } 
is uniformly in e bounded in L 2 ((0, T); H / 2 1 (f2)). Therefore it is possible to 
extract a subsequence of {e > 0} such that 

tf £ -»• weakly in L 2 ((0, T); W^fi)) 

as £ \ 0. 
Moreover, 

X^Dfo w £ ) - 0, x^/Vfl 6 - (4.1) 

as £ \ 0. 

Relabelling if necessary, we assume that the sequences converge them- 
selves. 
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On the strength of Nguetseng's theorem, there exist 1-periodic in y func- 
tions Q(x,t,y), P(x,t,y), U(x,t,y), Q(x,t,y), W(x,t,y), Q s (x,t,y) and 
U(x,t,y) such that the sequences {0 e }, {p £ }, {^ £ }, {l £ }, {w £ }, {V^ E } and 
{V x u £ } two-scale converge to Q(x,t, y), P(x,t,y), H(x,t,y), Q(x,t,y), 
W(x, t, y), V x -# + V y @ s (x, t, y) and V x u + V y U(x, t, y), respectively. 

Note that the sequence {div x w 6 } weakly converges to div x w and \u\ G 

o 

L 2 ((0, T); W\ (fi)). Last assertion for disconnected porous space follows from 

o 

inclusion d £ , \u £ \ G L 2 ((0,T); W\ (Q)) and for the connected porous space it 
follows from the Friedrichs-Poincare's inequality for u e and t? £ in the e-layer 
of the boundary S and from convergence of sequences {u £ } and {-# e } to u 
and $ respectively strongly in L 2 (Qt) and weakly in L 2 ((0,T); W^fi)). 

4.2. Micro- and macroscopic equations I.. 

Lemma 4.1. For all x G Q and y G Y weak and two-scale limits of the 
sequences {9 £ }, {p £ }, {^ £ }, {q £ }, { w£ }, {V x $ £ } and {V x u £ } satisfy the 
relations 



Q = ixq, Q = P + x(A)/e + 7 /); (4.2 

+ (1 - X ){div x u + div y U - ^(tf - (0> n ) + ls ) = 0; (4.3 

divyW = 0; (4.4 

w = x w + (i-x)u; (4-5 

e = x e + (i - (4.6 

q = p + (3of6 f + f\ (4.7 

^tt + (1 - m)(d«k« - ^(0 - (0) n ) + %) + (<%E/) Ys = 0; (4.8 

-J-tt + di^to - (1 - m)(^f (tf - (0) n ) - 7s ) + /? = 0, (4.9 



/? = ((div y U) Ys }a, f = (6)y f , 

1 

"77/ = (g) n - Pof(d J )n, (1 - m)7, = —(q)a - (3. 

Vo 

Proof. In order to prove first equation in (14.21) into Eq. fll.7p insert a test 
function i\f = erj) (x, t, x/e), where ip(x, t, y) is an arbitrary 1-periodic in y 
and finite on Yj function. Passing to the limit as e \ 0, we get 

V y Q(x,t,y) = Q, yeY f . (4.10) 



14 



The weak and two-scale limiting passage in Eq. (13.41) yield that Eq. (14.71) 
and the second equation in (14. 2p . 

Next, fulfilling the two-scale limiting passage in the equalities 

(i - x £ )p £ = o, (i - x £ )q £ = o 

we get 

{i-x)P = o, (i-x)Q = o, 

which justify first equation in (14.21) . 

Eqs. (14.31) . (14. 4p . (14. 8p . and (14.91) appear as the results of two-scale limiting 
passages in Eqs. (I3.2j) (13.41) with the proper test functions being involved. 
Thus, for example, Eq. (14.81) is just a subsequence of Eq. (l4.3p and Eq. (l4.9p is 
a result of two-scale convergence in the sum of Eq. (l3.3p and Eq. (l3.4p with 
the test functions independent of the "fast" variable y = x/e. Eq. (l4.4p is 
derived quite similarly if multiply the same sum of Eq. (13.31) and Eq. fl3.4p by 
an arbitrary function i\) £ = eif) (x, t, x/e) and pass to the limit as e \ 0. 

In order to prove Eqs. (14.51) and (14. 6p it is sufficient to consider the two- 
scale limiting relations in 



(1 - x £ )(w £ ~ u £ ) = 0, (1 - X £ W - $ £ ) = 0. 



Lemma 4.2. For all (x, t) G Qt an d y eY the relation 



□ 



<H{A (1 - X)(P(V, U) + D(x, u)) - (n + ^qx) • 1} = 0. (4.11) 
holds true. 

Proof. Substituting a test function of the form i\) e = eif) (x, t, x/e), where 
ip(x, t, y) is an arbitrary 1-periodic in y function vanishing on the boundary 
dfl, into Eq. (ll.7p and passing to the limit as e \ 0, we arrive at the desired 
microscopic relation on the cell Y . □ 

In the same way using additionally continuity equations (13.31) and (13.41) 
one gets from Eq. (11.81) 

Lemma 4.3. For all (x, t) G the relations 

A,0 S = 0, ye Y s: 



y i & s ' 

i: 

dn 

hold true. 



(4.12) 
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Now we pass to the macroscopic equations for the solid displacements. 

Lemma 4.4. Let p = mpf + (1 — m)p s , = (W)y>. Then functions 

u, , q, 7r, 9* , d satisfy in VLt the system of macroscopic equations 

T oPf-^r + T oPs{l-m)— - pF = (4.13) 
div x {X ((l - m)D(x, u) + (B(y, U)) Ys ) - (g + tt) • I}, 

def (3 2 d# fad* d(3 

r c p/ ^ + (r oCps + — )(1 - m)- - - - - fa- - (4.14) 

(1 - ™)A>»^ = x A{(1 - m)V x + (V y Q s ) Ys } + *. 

Proof. Eqs. fT4TT31 and fITOj) arise as the limit of Eqs.flEZD and (TOD with 
test functions being finite in an d independent of e. In Eq. (11.81) we have 
used continuity equations (13. 3p and (13 ,4p . □ 

4.3. Micro- and macroscopic equations II. 

Lemma 4.5. // pi = oo, then u = w and 9 = 

Proof. In order to verify, it is sufficient to consider the differences (u £ — w £ ) 
and (9 £ — $ £ ) and apply Friedrichs-Poincar'e's inequality, just like in the 
proof of Theorem 11.11 □ 

Lemma 4.6. Let p\ < oo and V = xdW/dt. Then 

dV 

T p f —-p f F = p l A y V-V y R-V x q, yeY f , (4.15) 

Wpf^ = x 1 /i 1 A 2/ 6 + ^^ + *, yeYf, (4.16) 
at m at 

du 

V=— , = 7?, ye 1 (4.17) 

for pi > 0, and 

dV 

T oPf^~ = ~ V y R - V x q + p f F } yeY f , (4.18) 

r c pf § = + *, y e y„ (4-19) 
at m at 

(xW -u)-n = 0, y e 7 (4.20) 

/or //1 = 0. 

In .Eg. (14.201) n is £/ie imit normal to 7. 
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Proof. Differential equations (14.151) and (14.181) follow as e \ from inte- 
gral equality (11.71) with the test function ip = Lp(xe~ l ) ■ h(x,t), where cp is 
solenoidal and finite in Yf vector-function. 

The same arguments apply for the Eq. (14.161) and Eq. fl4.19p The only one 
difference here is that we use the continuity equation (13. 3p to exclude the 
term x £ div x (dw e /dt). 

First boundary condition in (14.171) is the consequence of the two-scale 
i_ i 

convergence of {a^V x w £ } to the function \x\ V y W(x, t, y). On the strength 
of this convergence, the function W y W(x, t, y) is L 2 -integrable in Y. As 
above we apply the same argument to the second boundary condition in 
f T4TT7j) . The boundary conditions ffl~2T)|) follow from Eqs.flHD and < SM- ^ 

Lemma 4.7. If the porous space is disconnected, which is the case of isolated 
pores, then u = w. 

Proof. Indeed, in the case < jii < oo the systems of equations (14. 4p . 
( T4~T5|) and fHTTD . or (Ojl . (OBI) and fl4T20l) have the unique solution V = 
du/dt. □ 

4.4. Homogenized equations I. 

Lemma 4.8. If Hi = oo then w = u, 9 = d and the weak limits u, d, p, q, 
and 7r satisfy in the initial-boundary value problem 

ToP W + V(q + n)-pF= (421) 
div x {X A s : D(ar, u) + B s {div x u - + Bfq}, J 

(roc P + ^(l-m))^-^ + (at--)(%n = dw x {B 6 -W)+*, (4.22) 
Vo dt T] dt T] dt 

— 7T + C° : B(x, u) + a s (div x u - ^0) + a\q = 7, (4.23) 
Vo Vo 

-TT + dw x u + (1 ~ m)/3 °^ = 0, (4.24) 
Vo Vo 

q = p + ^QfTwd + mjf, (4-25) 

where the symmetric strictly positively defined constant fourth-rank tensor 
Aq, constant matrices Cq, Bq B{, strictly positively defined constant matrix 
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B and constants a , a{ and a s 2 are defined below by formulas (14.321) - (14.341) 
and (I4.37P and 

7 = K - -)(q)a ~ aS— Wn, -P=(l- m)—(0)n + -<?>n. 
Vo Vo Vo Vo 

Differential equations (I4.2ip and (14.221) are endowed with initial condi- 
tions at t = and x 6 Q 

(t + A) S )(tf - #o) = 0, r (u - u ) = r (— - v ) = 0; (4.26) 

and boundary conditions 

■&(x,t) = 0, u(x,t) = 0, x e S, t>0. (4.27) 

Proof. In the first place let us notice that u = w and 9 = •& due to Lemma 
1431 

The differential equations (14.211) follow from the macroscopic equations 
(I4.13p . after we insert in them the expression 

(D(y, U)) Ys = A{ : D(x, u) + B s (dW x u - ^(0 - (0) n )) + B{{q - (q) n ). 

In turn, this expression follows by virtue of solutions of Eqs. (l4.3p and 
(14. lip on the pattern cell Y s . Indeed, setting 

3 3 
U= J2 U i ^y)D ij + U (y)(dw x u - ^(0 - (0) n )) 

+U 1 (y)(q-{q)u) + U 2 (y){q) Q 



3 



n = a n ij (y) D a + My){&v x u - ^ - (tf) n )) 

+U 1 {y){q-{q) Q )+U 2 {y){q) n , 

1 dm duj 
ij 2 { d Xj + d Xi h 
we arrive at the following periodic-boundary value problems in Y: 

div„{(i - x )(P(v, u ij ) + J ij ) - if ■ i} = o, 

^ U ij + (i _ x )div y UV = 0; 



where 



(4.28) 



div y {A (i - xWv, u ) - n Q • 1} = 0, 

i ; n + (l- X )(div J/ t/ + l)=0; 

1 



no 



dw y {X (l - xWv, Ui) - (Hi + -x) ■ 1} = o 
-i-IIx + (1 - X )^yUx) = 0. 



div,{A (i - u 2 ) - (n 2 + -x) • 1} = o, 

1-yK,.. . 1 



1 (1 - y) 

Vo (1 - m) 



«div„E7 2 )y. + 



0. 



Note, that 



(4.29) 
(4.30) 

(4.31) 



/3 = V(div^>y s <A^ + (div,C/ )y s (div^ - - (0> n )n+ 

(div„aT)y,(g - (g) n ) n + (div y J7 2 )y s (g) n = (div y E/ 2 )y s (g) n 

due to homogeneous boundary conditions for u(x,t). 

On the strength of the assumptions on the geometry of the pattern "liq- 
uid" cell Y s , problems (I4.28p ~ (14.311) have unique solution, up to an arbitrary 
constant vector. In order to discard the arbitrary constant vectors we de- 
mand 

(U ij ) Ys = (U ) Ya = (Ufa = (U 2 )y s = 0. 

Thus 

3 3 

k = E JlJ ® Jl3 + A i> A i = E <(* - ® Jii - ( 4 - 32 ) 

Symmetry and strict positiveness of the tensor Aq have been proved in [TO] . 

Finally, Eqs. (147231 - (gi25} for the pressures follow from Eqs. (jiT7fl - (Ojl . 
after we insert in them the expression 

(div^)y =C S :D(x,u) + ^(div :cW -^(^--W c )) + a s 1 (g-(g) n )+^(g) n 
where 

3 

I?* = (D(y, t/ ))y , 5? = (0(y, E^y,, C s = ]T (div^y.^, (4.33) 
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a s = (dw y U } Ys = a s -l+m, a{ = (div y Ui)y s , a* 2 = (dw y U 2 } Ys - (4.34) 
Now for i = 1,2,3 we consider the model problems 

A y ef = 0, ye Y s , (4.35) 

— — = -e t n, y G 7 

on 

and put 

3 

1=1 

Then B s solves the problem (l4.12l) -(??) and if we insert an expression (V y O s )y s 
into fl4TT4"| we get 

3 

B e = x 0s ((l - m)l + ^(V y e^)y s ® a). (4.37) 

i=i 

All properties of the matrix B 9 are well known ( see [2], [7]). □ 

Lemma 4.9. // the porous space is disconnected, then w = u and the weak 
limits 9? , u, i? ; p, q, and n satisfy in Qt equations ( I4.2ip . (I4.23[) . (14.241) . 
(I4.10p . where Aq, Cq,Bq B{, B^ , a s , a\ and a s 2 are the same as in Lemma 
\4-8[ the state equation (14 .7p . and heat equation 

80* (31 M (3 0s dn 1 dq 

roC pf ^ + (roc ps + -)(1 - m)- - — - + (a, - _)<_>„ = 

div x (B e ■ V#) + 

(4.38) 

where for ^i > and r > 

0/(:M) = mtf(cM) + fb 6 {t-r){— (^L^- + ^) -^)(x,r)dr. (4.39) 
If Hi > and t = 0, then 

9 f (x, t) = m#(x, t) - c e A^^-(t) + *)). (4.40) 
Finally, if fii = 0, taen 

^(aj,t) = mi? (aj) + — f\—^(r) + ^(x,r))dr. (4.41) 

i/ere and are defined below by formulas (I4.42D — (14.44j> . 

The problem is endowed with initial and boundary conditions (14.261) (14.271) . 
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Proof. The only one difference here with the previous lemma is the heat 
equation for $ and the state equation for pressures, because 9 ^ The 
function 9* = (@)y / now is defined from microscopic equations (I4.16P and 
(I4.17p . if fMi > and microscopic equations (I4.19p . if H\ — 0. 

Indeed, the solutions of above mentioned problems are given by formulas 

G = #(x,t)+ [ e{(y,t-r)h(x,r)dT, 
Jo 

if yUi > and r > and 

m at 

if \i\ > and r = 0, where 

and functions 0[ and Bp are 1-periodic in y solutions of the problems 

T ° c pf~df = *m A y®i, V^ Y h I (4.42) 

9{(y,0) = l, yeY f ; 6{ = 0, ye 7, J 

and 

xi^ 1 A 1/ eJ = l, t/6 7/; 6^ = 0, 2/G7- (4.43) 

Then, in accordance with definition, the function 9* is given by (14.391) or 
fr4T4T?|) . where 

b]{t) = (e{) Yf , c] = (e %. (4.44) 

If /ii = 0, then is found by a simple integration in time. □ 
4.5. Homogenized equations II. 

Let fi\ < oo. In the same manner as above, we verify that the weak limit 
u of the sequence {u 6 } satisfies some initial-boundary value problem likes 
problem (14.211) - (I4.27P because, in general, the weak limit w of the sequence 
{w 6 } differs from u. More precisely, the following statement is true. 

Lemma 4.10. // [i\ < oo then the weak limits u, u>f , 9*, d, p, q, and n of 

the sequences {u £ }, {x £ w £ }, {x £ 9 e }, {$ £ }, {p £ }, {q 6 }, and {n £ } satisfy the 
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initial-boundary value problem in Qt, consisting of the balance of momentum 
equation 

To(Pf-7^ + Ps(l-m)—) + V(q + ir)-pF= ^ 
div x {\ AQ : D(x, u) + B^diVxU + Bfq}, 

where Ag, Bq and Bf are the same as in (I4.2ip . continuity equation (14.231) . 
continuity equation 

-(tt + (g)n) + div x w f + (1 ~ m)/3 ° s (^ - (#) n ) = (m - l)div x u, (4.46) 

?7o ^0 



siate equation (14.71) . aeat equation (I4.38j) and Darcy's law in the form 
dw? du /"* „ / . , _ „ S 2 1 

i/r > and /ii > 0, Darcy's law in the form 
dw? du 



du /"* <9 2 w 
— + / S 1 (// 1 ,t-r)-(-V 3: g + p / F-rop / -^-)(a ! ,r)dT (4.47) 



S 2 (mi) • (-V^ + p/F) (4.48) 



^ dt 

if T o — and, finally, Darcy's law in the form 

r) w f f) u 1 /•* 

= B 3 - — + (ml-B 3 ). / (-V«g + p,F)(*,r)dT (4.49) 

at at r pf Jo 

if Pi = 0. The problem is supplemented by boundary and initial conditions 
(I4.26P - (14.271) for the displacement u and temperature $ of the rigid compo- 
nent and by the boundary condition 

w f (x,t) ■ n(x) = 0, (x,t) E S = dtt, t>0, (4.50) 

for the displacement of the liquid component. In Eqs. (14.45)) - (14.491) n(x) 
is the unit normal vector to S at a point x e S, and matrices Bi(p 1 ,t), 
S 2 (/ii) ; and B% are defined below by formulas (14.511) - (14. 56)) . 

Proof. Eqs. (I4.45P and (I4.46P derived in a usual way like Eqs. (14.211) and 
(14.24)) . For example, to get Eq. (l4.46p we just expressed div x w in Eq. (14.91) 
using homogenization in Eq. (14.51) : w = + {l—m)u. Therefore we omit the 
relevant proofs now and focus only on derivation of homogenized equations 
for the velocity v in the form of Darcy's laws. The derivation of Eq. (14.501) 
is standard 
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) If Hi > and r > 0, then the solution of the microscopic equations 
(B3D, QHE>} and (l4TTTj) is given by formula 

where 

3 

i=i 

and functions V l (y,t) are defined by virtue of the periodic initial-boundary 
value problem 

roP/^-ZiiAV' + VQ^O, div tf V* = 0, yeY f ,t>0,\ (451) 
V l = 0, yG 7 ,t>0; r oPf V i (y,0) = e l , y eY f . J 

In (I4.5ip e« is the standard Cartesian basis vector of the coordinate axis £j. 
Therefore 

S 1 (/x 1 ,t) = (B{( 1 /,t))y.. (4.52) 

b) If To = and fix > then the solution of the stationary microscopic 
equations ( 14.41) . ( 14.151) and ( 14.171) is given by formula 

V=— + B f 2 (y)-(-Vq + p f F), 

where 

3 

B f 2 (y) = J2u\y)®e l , 

i=l 

and functions U l (y) are defined from the periodic boundary value problem 



-fiiATP + Vi? = e i; divyET* = 0, ye Y f , 
U* = 0, ye!. 



(4.53) 



Thus 

BM = (B(((y))y 3 . (4.54) 

Matrices Bi(fii,t) and ^(/ii) are symmetric and positively defined [T 
Chap. 8]. 
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c) Finally, if r > and fx\ = 0, then in the process of solving the system 
(14.40 . 04.181) and 04. 201) we firstly find the pressure R(x,t,y) by virtue of 
solving the Neumann problem for Laplace's equation in Yf. If 

h(x, t) = -r p f — (x, t) - Vq(x, t) + p f F(x, t), 

then 

3 

R(x, t, y) = ^ Ri(y)ei ® h(x, t), 

i=l 

where R l {y) is the solution of the problem 

ARi = 0, ye Y f ; VRi n = ne u y e 7. (4.55) 

The formula 04.491) appears as the result of integration with respect to time 
in the homogenization of Eq. 04.181) and 

3 

B 3 = J2(VRi(y))Y s ® e h (4.56) 

i=i 

where the matrix (ml — ^3) is symmetric and positively definite [T4"l Chap. 
8]. □ 
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